Comment on "Negative Heat Capacity for a Cluster of 147 Sodium Atoms" 
by Schmidt et al. 

In a recent letter [I], Schmidt et al. analyze a region where the specific heat of a Naf 47 
cluster becomes negative. In Boltzmann-Gibbs (BG) statistical mechanics the specific heat 
can never be negative in the canonical ensemble. However, as the authors are well aware, no 
such property is valid in general for the microcanonical ensemble. This is indeed the case of a 
variety of physical situations ranging from self-gravitating systems to melting atomic clusters 
and fragmenting nuclei (see [p]] and references therein). An important feature that all such 
systems have in common is thermodynamic nonextensivity, reflected, as also mentioned 
in in the fact that the energy associated with N particles is not proportional to N. 
Nonextensivity occurs everytime the range of the interactions is not negligible compared to 
the (linear) size of the system. Two important physical examples are: 

- short-range interactions in small systems 

- long-range interactions in systems of any size [§-[§• 

Now, if the system is isolated, nonextensive, in an equilibrium-like state, and exhibits, for 
some energy range, a negative specific heat, then not one but at least two possibilities must 
be considered: 

1) the system is in its thermodynamically stable state and then the BG microcanonical 
ensemble is to be used, as preconised in || and done in 0. 

2) the system is in some kind of dynamic metastable state, and then some other thermo- 
statistical approach might well be necessary. One such alternative description is nonextensive 
statistical mechanics 0, which has proved to be useful in a variety of complex situations 
(e.g., fully developed turbulence ||, electron-positron annihilation producing hadronic jets 
0, motion of Hydra viridissima ||). 

The metaequilibrium state to which we make reference as a second possibility has indeed 
been observed in magnetic-like Hamiltonians such as arrays of long-range coupled planar 
rotators (the so called HMF || and its generalization, the a — XY model [TCfl), among 
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others, including fluid-like systems [[LI]]. To illustrate this point, we report in fig. 1 a 
numerical simulation typical of the relaxation process in the HMF model. As a function 
of time, two plateaux are observed for the average kinetic energy per particle and the BG 
formalism is violated for the first plateau (panel a), where we find a region of negative 
specific heat (panel c) ]Tj| . More precisely, the velocity distribution is well described by the 



standard Maxwellian in the late plateau (panel b), whereas for the early plateau it appears to 
be described by the type of distribution emerging within nonextensive statistical mechanics 
18| . Similar results and anomalies have been recently described in other models, including 



first order phase transitions [|l2Hl4l 



The relevance of long range interactions in sodium clusters and similar atomic liquids is 
well known [[nj. Consequently metastable states could very well occur in relatively small 



clusters such as those studied in |j] , as in various colloidal glasses [|TIJ and other clusterized 



systems |[L7|| . For example, in colloidal glasses it is common belief that the structural arrest 
at the glass transition is due to the increasingly size of strongly cooperating sets of atoms 
| |16| |. Under these circumstances, and without further experimental evidence, application of 
BG equilibrium statistical mechanics in might have no support. Since the analysis of 
the authors only applies to one among at least two physically important possibilities, the 
theoretical interpretation they draw from their experiments JTJ is not necessarily true. 
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(Brazilian agencies) is also acknowledged. 

C. Tsallis 1 ' 2 , B.J.C. Cabral 1 , A. Rapisarda 3 and V. Latora 3 

1 Centro de Fisica da Materia Condensada, Universidade de Lisboa,Av. Prof. Gama Pinto 
2, P-1649-003 Lisboa, Portugal 

2 Centro Brasileiro de Pesquisas Fisicas, Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro- 
RJ, Brazil 

3 Dipartimento di Fisica e Astronomia, Universita di Catania, and INFN, Sezione di Catania, 



Corso Italia 57, 1-95129 Catania, Italy 

* Email: tsallis@cbpf.br, ben@adonis.cii.fc.ul.pt, andrea.rapisarda@ct.infn.it, 
vito.latora@ct.infn.it 

PACS numbers: 05.70.Fh, 36.40.Ei 



3 



REFERENCES 

[1] M. Schmidt, R. Kusche, T. Hippler, J. Donges, W. Kronmuller, B. von Issendorff and 
H. Haberland, Phys. Rev. Lett. 86, 1191 (2001). 

[2] M.E. Fisher and J.L. Lebowitz, Commun. Math. Phys. 19, 251 (1970). 

[3] C. Tsallis, J. Stat. Phys. 52, 479 (1988); S. Abe and Y. Okamoto (eds.), Nonextensive 
Statistical Mechanics and Its Applications, Series Lecture Notes in Physics (Springer- 
Verlag, Heidelberg, 2001); For a regularly updated bibliography of the subject see 
http:/ /tsallis. cat. cbpf.br/biblio.htm . 



[4] J. Barre, D. Mukamel, S. Ruffo, Phys. Rev. Lett. 87,030601 (2001). 

[5] D.H.E. Gross, Microcanonical thermodynamics: Phase transitions in "small" systems, 
Lecture Notes in Physics 66 (World Scientific, Singapore, 2001). 

[6] C. Beck, Phys. Rev. Lett. 87, 180601 (2001) T. Arimitsu and N. Arimitsu, Phys. Rev. 
E 61, 3237 (2000). 

[7] I. Bediaga, E.M.F. Curado and J. Miranda, Physica A 286, 156 (2000); C. Beck, Physica 
A 286, 164 (2000). 

[8] A. Upadhyaya, J.-P. Rieu, J. A. Glazier and Y. Sawada, Physica A 293, 549 (2001). 

[9] V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett. 80, 692 (1998); Physica D 131, 
38 (1999); Phys. Rev. Lett. 83, 2104 (1999); V. Latora and A. Rapisarda, Nucl. Phys. 
A681, 406c (2001). 

[10] C. Anteneodo and C. Tsallis, Phys. Rev. Lett. 80, 5313 (1998); A. Campa, A. Giansanti, 
D. Moroni and C. Tsallis, Phys. Lett. A 286, 251 (2001). 



[11] E.P Borges and C. Tsallis, Physica A (2002) in press ||cond-mat/ 0109 504 



[12] Y. Sota, O. Iguchi, M. Morikawa, T. Tatekawa and K.I. Maeda, Phys. Rev. E 64, 056133 
(2001). 



[13] Lj. Milanovic, H.A. Posch and W. Thirring, Phys. Rev. E 57, 2763 (1998). 
[14] A. Torcini and M. Antoni Phys. Rev. E 59, 2746 (1999). 
[15] J.P.K. Doye and D.J. Wales, Science 271, 484 (1986). 

[16] E.R. Weeks, J.C. Crocker, A.C. Levitt, A. Schofield and D.A. Weitz, Nature 287, 627 
(2000); see also G. Foofi et al. [|cond-mat/0111221| . 



[17] B.J.C. Cabral and J.L. Martins, J. Non- Crystalline Solids (2002), in press. 



18] V. Latora, A. Rapisarda and C. Tsallis, Phys. Rev. E 64, 056134 (2001) and ||cond- 
mat/0109056|| ; A. Campa, A. Giansanti and D. Moroni, ||cond-mat / 1 09 1 TSfl . 



5 



FIGURES 




Time 



6 



FIG. 1. Time evolution of twice the average kinetic energy per particle ("temperature") for a 
HMF system with N=500 (red circles) at energy per particle U = 0.69. The corresponding velocity 
distribution is not Maxwellian for the first, metastable plateau (panel (a)), whereas standard 
statistical mechanics is valid for the second, stable plateau (panel (b)), open circles. The Gaussian 
canonical prediction is also shown as a full curve. Correspondingly, we find a negative specific 
heat, Cy = dU/dT, in the metastable state (Quasi-Stationary State (QSS) regime) , as shown in 
panel (c) (open circles) in comparison with the canonical caloric curve, while the specific heat for 
this hamiltonian is always positive when BG equilibrium is reached (panel (d)). Such anomalies 



are strongly related to nonextensivity and have been found in various systems |ll|,|l2,14,li 
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